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^ Abstract 

^ It is consistent with the measurement of ^13 ~ 0.15 by Daya Bay to suppose that, 

*""' in addition to being unitary, the neutrino mixing matrix is also almost hermitian, and 

^^ thereby only a small perturbation from diag(+l, —1, —1) in a suitable basis. We suggest 

this possibility simply as an easily falsifiable ansatz, as well as to offer a potentially 



(N 



^ 
O^ 



useful means of organizing the experimental data. We explore the phenomenological 
implications of this ansatz and parametrize one type of deviation from the leading order 
relation \Ve2\ ~ |^i|- We also emphasize the group-invariant angle between orthogonal 



Q matrices as a means of comparing to data. The discussion is purely phenomenological, 

,-S^ without any attempt to derive the condition yT ~ y from a fundamental theory. 

cn 
> I A Phenomenological Ansatz 

f — 

\^ The neutrino mixing matrix V is defined by v^ = VaiJ^i, where a = e,iJ,,T denotes the 
CN charged lepton mass basis ("flavor basis"), and i = 1,2,3 denotes the neutrino mass basis. 

cn 

o 



en The relevant part of the Lagrangian written in the flavor basis reads 
O 



£ = - ^ rriaeaea - \va{Mu)apVp + h.c. (I.l) 



c3 the nii real and positive. 

Assuming that the three light neutrinos of the Standard Model are Majorana, the mag- 
nitudes of the entries of V are constrained by data to be 

/0.78-0.84 0.52-0.61 0.13-0.17\ 
iKxpl ~ 0.40-0.58 0.39-0.65 0.57-0.80 . (1.2) 

\0.19-0.43 0.53-0.74 0.59-0.81/ 



To obtain Eq. (1.2) we have used 0.550 < ^12 < 0.658 and 0.620 < ^23 < 0.934 from the 



work of Gonzalez-Garcia, Maltoni, and Salvado pQ, and 0.135 < ^13 < 0.171 from the recent 

1 



results of Daya Bay [2J. The ranges in Eq. (1.2) are correlated in such a way as to preserve 
the unitarity condition V'^V = I. 

In an effort to obtain a theoretical understanding of the mixing matrix, one might sup- 



pose that the numerical values in Eq. (1.2) arise as a small perturbation from a "simple" 
ansatz. As a straw man argument for what such an ansatz might be, consider an older global 
best fit given by [3] 
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Kxp.oidl- 0.22-0.56 0.44-0.73 0.57-0.80 . (1.3) 



Simply by looking at the ranges in Eq. (1.3), we observe that it was once numerically consis- 
tent to suppose that V is hermitian. Without any theoretical motivation, we now suppose 
that the true mixing matrix satisfies the leading order relation V'' ~ V, and then we study 
small deviations required to fit the new data. We propose this rather ad hoc constraint in 
the spirit of trying to make sense of the data by reducing the number of free parameters in 
the neutrino sector [4J. This exercise is intended partially to illustrate that there are still 
many possibilities for what the true mixing matrix might be. 

II Real Symmetric Mixing Matrix 

As a warmup, we first consider the case for which V is real. Then V is orthogonal, meaning 
V'^V = I, and our hermiticity ansatz amounts to imposing the condition V'^ = V. The old 
experimental bounds adjusted for compatibility with this ansatz are 



'0.77-0.86 0.50-0.56 0.21-0.22^ 

exp, old I 



IKxTomI - I 0.50-0.56 0.44-0.73 0.57-0.71 I . (II.l) 



,0.21-0.22 0.57-0.71 0.59-0.82, 



Again, the ranges are correlated, as required by V'^V = I [See Eq. (II.4)]. Immediately we 
see that the Daya Bay observation of \Ve3\ < 0.21 rules out this ansatz as an exact predic- 
tionj but otherwise it is still consistent with Eq. (1.2). 



If V is real symmetric, then it can be diagonalized by an orthogonal transformation: V = 
XdX^, where d is diagonal and X is orthogonal. Then V'^V = I implies V"^ = Xd'^X^ = I, 
so that d"^ = I. Thus our ansatz amounts to proposing that, in a particular basis, the neu- 
trino mixing matrix is diagonal with entries equal to ±1. 

We now have a choice as to how to arrange the minus signs in d. Two of the nonzero 
entries in d must have the same sign, while the third must have a sign opposite to that of the 



-'^Much of this work was completed before the measurement of nonzero reactor angle, 6*13 ^ 0.15. The fact 
that the reactor angle is relatively large, meaning closer to ~ 0.2 than to zero, is what keeps our analysis 
relevant. 



first twcQ In other words, we get to choose which 2-diniensionaI subspace of d is proportional 
to the identity matrix. This choice is arbitrar}Q to fix the discussion, we choose 



ci= -1 (II.2) 




so that d equals minus the identity matrix in the (2, 3)-plane. 

Any rotation matrix in 3 dimensions can be written as a product of independent rotations 
about each of 3 mutually orthogonal axes. That is, given the rotation matrices 

/I \ / C2 S2\ / Cs ^3 0\ 

X,= \0 Ci -Si] , X2=\ 1 , X3 = -53 C3 (11.3) 

\0 Si Ci I \-S2 C2J \ 1/ 

where Cj = cos (pi and Si = sin if j, we can write X as a product of the three Xj in any 
ordeiQ Since d is proportional to the identity matrix in the (2, 3)-plane, it is unchanged by a 
rotation about the first axis: XidX^ = d. Thus one of the parameters in V^ = XdX'^ drops 
out, leaving us with a two-parameter ansatz for the mixing matrix. Choosing the ordering 
X = X3X2X1 implies 



C|cos(2¥;3)- 


-SI 


-C|sin(2^3) 


— C3 sm{2ip2 


X 




-CI cos(2^3) - SI 


S3 sm{2ip2) 


X 




X 


— COs(2y92) 



V = \ X -CI cos(2^3) - ^2' ^3 sin(2^2) (11.4) 

\ X X — COs(2v92) / 

where we have displayed only the upper triangle of V since by construction it is symmetric. 



The values for the angles consistent with Eq. (1.3) turn out to be 0.32 < v^2 < 0.42 and 



1.20 < V93 < 1.27, where the angles are expressed in radians, as shown in Fig[T} 

As an arbitrarily chosen "typical" example of hermitian mixing, the values {ip2, 'Ps) 
(0.35, 1.23) imply a mixing matrix 




(11.5) 



^The solution d ~ I trivially satisfies cP = /. This would result in y = XdX^ = XX^ — I, which is of 
course experimentally unacceptable. 

•^For example, let X — X'P, where X' is orthogonal and P is a permutation matrix. Then V = XdX^ = 
X' d'{X'Y^ is of the same form as before, but with a new diagonal matrix d' — PdP^ with the two minus 
signs permuted. Of course, we are also free to multiply V and hence d by an overall sign. 

'^At this point we should emphasize that ipi are not the three PMNS angles that parametrize the mixing 
matrix V — XdX^ . That is why we have chosen to denote their sines and cosines by capital letters, in 



contrast to the notation in Section III for the usual PMNS angles. 
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Figure 1: The values of ip2 and (p^ in the quadrant ((/?2, Vs) ^ ([Oj f ]> [Oj f ]) consistent 



with an older set of oscillation data, Eq. (1.3). This may be a useful starting point 



about which to perturb in order to fit the new data, Eq. (1.2) 



where we have rearranged the minus signs into a standard form[j Compare this with the 
often-studied "tribimaximal mixing" ansatz [HI E] 
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(11.6) 



These two matrices appear quahtatively "very different," given that one has K3 ^0.22 while 
the other has Ve3 = 0. To make this notion more precise, definajthe 5'0(3)-invariant angle 
between two special orthogonal matrices V and V: 



Q{yy) = cos-Wii{y^v')) 

The matrices Vo.35,1.23 and Vtb are separated by an angle 



0(K).35,1.25,Vtb) ~ 



0.20 ~ 11° 



(11.7) 



(11.8) 



in S'0(3). As a related example, one might compare to another ansatz with the same atmo- 
spheric and reactor angles as tribimaximal mixing (^23 = f and 6*13 = 0, respectively), but 



III 



^So as not to interrupt the logical flow we will postpone discussion of rephasing V until Section 
^Another measure ^2] of the 5'0(3)-invariant distance between matrices is 'D{V,V') — itr(/ — V V'). 

Here we choose the angular distance because it provides an intuitive notion of "large" versus "small" in terms 

of an angle 8 ranging from to 7r/2. 



with the solar angle related td^ the golden ratio: 6*12 
[3 El [91 [ini HI] The PMNS matrix for this ansatz is 



tan"^(l/v9), with ip = |(1 + V5 ). 



K 



golden 




0.53 
0.60 0.71 
0.60 -0.71, 



(11.9) 



This is separated from the matrix of Eq. (II. 5) by an angle 6(Vo.35,i.25, V^oWen) ~ 0.20 ~ 11° 



approximately the same as for Vtb- More generally, we see that the entire family of "/ir- 
symmetric" mixing matrices is approximately separated from the entire family of hermitian 
mixing matrices by ~ 11° in 5*0(3). 



We now know, due to Eq. (1.2), that neither Eq. (II. 5) nor Eq. (II. 6) is correct at low 



energy, but both may still serve as leading order predictions. The present oscillation data 
therefore admit two possible starting points that are separated by an angle ~ 11° in the set of 
all possible 3-by-3 special orthogonal matrices. This is simply intended to give a quantitative 
measure of uncertainty in our theoretical understanding of the mixing matrix. 



We might also like to obtain a quantitative sense of how different the new data is from the 
old. As perhaps an overly simplistic approach [12], we extract the arithmetic mean values 



for the matrices in Eqs. (1.2) and (1.3) and fix the signs according to a chosen convention: 



T^mean 
exp 



-0.810 
0.490 
0.310 



0.565 
0.520 
0.635 



0.150 

0.685 

-0.700, 



T^mean 
exp, old 



-0.815 


0.565 


0.110 


0.390 


0.585 


0.685 


0.380 


0.555 


-0.705 



(11.10) 



Neither of these matrices is orthogonal. To correct for this, we define "corrected" versions 
of these matrices by multiplying on the right by some other matrix F yet to be determined: 



— TAmeanp Thpn ( V/'t^'^^^^T 'W'mca.n 



•^exp ' exp -" ■ -" '■'-^^■^ \ • Qxp / ' exp 

real symmetric, so it can be written as (V^ 



T^iy~YV~V. The matrix {y^^^^ 






IS 
exp / ' exp -" ■ ■" '■'■"^ ij-i^v-j-ij^ y^ . gjjp J • gjjp iu 

™^^'^^^1/?r'' = Sk^S^, where A^ is diagonal with 



exp / exp 



real positive entries, and S'^ S = I. Thus if we fix^ 



r 



SA-^S^, we have {V~fV~ 



I. 



The "corrected" experimental mixing matrix 



T^^mean -- T^mean C A ^1 C^ 
*^exp *^exp 'J-'*- "J 



is orthogonal. Carrying out this procedure for V^^p*^" and Kxpfoid gives 



(11.11) 



T^inean 
exp 



-0.811 0.565 
0.494 0.525 
0.312 0.637 



0.151 

0.693 

-0.705, 



K 



exp, old 



-0.824 0.564 0.053 
0.437 0.574 0.693 
0.361 0.594 -0.719, 



(11.12) 



^The references discuss various different proposals for relating the solar angle to the golden ratio. We 



arbitrarily choose the particular implementation of Eq. (II. 9 1 to be concrete. 

®The ordering of the three eigenvalues in A^ and the corresponding eigenvectors in S is arbitrary and does 
not change the m atrix T = SA~^S'^ . This is conceptually the same freedom as that of arranging the signs 



in d of Eq. (II. 2 1 and the columns in X without changing V = XdX 



Using these, we find 

0(V;~K^roid)^ 0.08-5° (11.13) 

Thus, in some quantitative sense, the old data is "close" to the new data and therefore may 
help characterize possible ansatze away from which the true mixing matrix might be only a 
small perturbation. 

Ill Complex Hermitian Mixing Matrix 

We are now ready to consider a fully complex mixing matrix V. If V is hermitian, then it 
can be diagonalized by a unitary transformation: V = XdX\ where d is diagonal and X is 
unitary. The unitarity condition V'^V = I and the hermiticity ansatz V"^ = V imply d"^ = I, 
just as for the real case. Again we take d = diag(+l, —1, —1). 

The most general unitary 3-by-3 matrix has 3^ = 9 independent real parameters. We 
now briefiy recapitulate the justification behind the standard angular parameterizatioE|j of 
a unitary matrix [131 CI]- If t^ is a unitary matrix, then a matrix V whose elements are 
Vij = e^^^^'^y^^'Uij is also unitary. The angles Xi and yj together constitute 5 independent 
parameters, not 6, since they enter only in the combination x, + i/j. Three of the four re- 
maining parameters in U can be taken as the three independent rotations from Section |TT} 
The final parameter can be included as a non-removable phase in one of the rotation matrices. 

Thus we arrive at the usual PMNS parameterization of the neutrino mixing matrix, V = 
ICVpumM, where /C = diag(e^'=% e^'^^e^'^^)^ _A/f ^ diag(e*'',e^'", 1) and 

sise-^-^cpX /-ci2 si2 0\ 

VpMNS = I C23 S23 I I 1 | ( S12 C12 (III.l) 

1/ 

where cjj = cosOu and su = sinO ij. The matrix /C is unphysical and can be chosen 
arbitrarily. The "Majorana" matrix Ai is physical if the neutrinos are Majorana, but it 
drops out of oscillation probabilities and hence is not observable in oscillation experiments. 
The four parameters 6jj and 6cp contribute to oscillations. With the chosen sign conventions, 
the PMNS matrix is 

-C13C12 C13S12 Si3 

^PMNS = I C23S12 + S23S13C12 C23C12 - S23S13S12 S23C13 I , Si3 = Sige^'^^P (III.2) 

,'523'5l2 ~ C23S13C12 S23C12 + C23S13S12 — C23C13 , 





which explains the signs in Eqs. (II. 5) and (II. 6). 



As discussed in Section |TT| we are free to choose the sequence in which the three rotations 



in Eq. (II.3) occur. Moreover, we are free to insert the phase angle 6cp into any of the three 



^Different parametrizations of the neutrino mixing matrix have been studied in the literature [1511161 ITT] . 



rotation matrices in Eq. (II. 3). The PMNS parameterization (III.l III. 2) is useful for neu- 



trino mixing since the relation Ve3 = S13 e ^^'^^ makes clear that 6cp drops out of the mixing 
matrix if Ves = 0. 



Even if the above parameterization is convenient for V, we should determine the param- 
eterization most convenient for the unitary matrix X, which we remind the reader is defined 
by 1/ = XdXl 

Since we chose d = diag(-|-l, —1, —1) to be proportional to the identity matrix in the 
(2, 3)-plane, we should put the analog of ^cp into a rotation about the first axis. We write 
X = N*XsX2XiK, where A^ = diag(-e^^ +e^^ +1), K = diiig{e'''\e'''\e'''^), X2 and X3 



are given in Eq. (II. 3), and finally 





Xi = I Ci -Si e-'"^ 

,0 516+*" Ci 



(III.3) 



Then both K and Xi drop out of the mixing matrix V = XdX'^ . Putting this into the form 
V = ICVpMNsM give(3 



V, 



PMNS 



^Cl cos(2v93) - SI Cl sin(2^3) C^ sin(2^2) 

-Cl cos(2^3) - Si sin(2(^2)53 



X 
X 



(III.4) 



X — COs(2(/92) 

which is precisely the same PMNS matrix as for the real symmetric mixing matrix ansatz 



in Eq. (II.4) up to rearranging the signs. When verifying that the signs are correct, keep in 



mind that sin(2(y93) > but cos(2(y93) < 0, and | cos(2(/93)| > sin(2(/93) for the allowed values 
of (/33 (see Fig. [l|)- 

Thus the ansatz of hermiticity V'^ = V, if it were exactly true, would predict that neu- 
trino oscillations conserve CP. 



IV Deviations from Hermiticity 

As can be seen from the data in Eq. (|I.3l), the measurement of |VJ=3| < 0.21 has immedi- 



ately ruled out the possibility that V is exactly hermitian. Moreover, as shown in Eq. (III.4), 
measuring any CP violation in neutrino oscillations would also signal deviations from V^ ~ V. 



In either case, V could still be approximately hermitian up to small corrections, so we should 
find a way to write deviations from hermiticity as a perturbation expansion in a small param- 
eter. In any such parameterization away from V"^ = V, we obviously still need to maintain 



^"More explicitly, the matrix N can be written as A^ = ~dA4, where d is the matrix of Eq 
and Ai is the matrix of Majorana phases defined above Eq, 



V = XdX'f = M*idX3X2dX^{Xld)M is put into the form V 
and /C = M*. 



{11.2) 



(III.l), so that r — p and s ~ a. Then 
ICVpMNsM with VpMNS = dX:iX2dX2 X'^ d 



the unitarity condition V^V = I. 

One possibility is to consider deviations from d'^ = d while maintaining d'^d = I. Consider 
the orthogonal matrix V^ = {X3X2)de{X3X2)^ , where 

/I 

(ig = —cose H-sine . (IV.l) 



In the limit e — ?■ 0, this matrix recovers the most general real, orthogonal, symmetric 3-by-3 
matrix given in Eqs. ( |II.4[ ) and ( |III.4| ). We have V = V{e = 0) + e6V + O(e^), where[^ 



6V 






-S2 


+C2S3 


4-^2 





+C2C3 


C2SZ 


—C2C3 






(IV.2) 



manifestly parameterizes a particular type of deviation from the real symmetric ansatz. In 
particular, it parameterizes a decrease in |Ve3| and a corresponding increase in |K-i|- 



Now we are ready to confront the new data of Eq. (1.2). The hermitian ansatz (III. 4) is 
not compatible with |Kxp| due to the tightened bounds on V^a. As a concrete example, recall 



the unperturbed matrix of (II. 5) with (y92, V's) = (0.35, 1.23), which we repeat for convenience: 

/-0.80 0.56 0.22 \ 
^o.35,i.23(^ = 0) ^ 0.56 0.57 0.61 . (IV.3) 

\ 0.22 0.61 -0.76/ 

Perturbing this matrix in the form V = V{e = 0) + e6V with e = 0.074 and 6V given in Eq. 



(IV.2) implies a mixing matrix 



-0.80 0.58 0.15 
V^ { 0.53 0.57 0.63 
0.28 0.58 -0.76. 



(IV.4) 



which is compatible with the bounds in |V"exp|. 

Since we are breaking the hermiticity ansatz, we re-introduce the possibility of CP viola- 



tion. Thus we can generalize Eq. (IV.l) to the case 



de,5 



'1 

— COS e 
, — e~*'^sin£: 



+e 





+15 



sme 



(IV.5) 



cose 



which implies a CP angle 



3CP 



tan 



& sin 6 sin e 



'I + cose)C3S2 + S'scos^sine 



tan(y?3) sin^ . 

2^2 ^ 



(IV.6) 
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Here we find it convenient to keep the signs as in Eq. (II.4I, and then to convert to the signs in Eq 



(III.4) only after adding V{e = 0) + e6V. 



For the previous case ((y92,v23,£) = (0.35,1.23,0.074) this gives (5cp ~ 0.30 sin 5. The basis- 
independent Jarlskog invariant J = lm.{VeiV*2V*iV^2) is given to leading order in e by 

J= -lC2sin(4(^2)sin(4(^3)£ sin5 + 0{e^) (IV. 7) 

where the coefficient of e sin 6 ranges from 0.217 to 0.229 over the range of values for {ip2, V^s) 
given in Fig. [Tj 

The purpose here is simply to illustrate that although present data falsify the ansatz V'^ = V, 
it is entirely possible that the true mixing matrix is only a small perturbation away from 
being hermitian. Just as the CKM matrix is a small perturbation from the identity, the neu- 
trino mixing matrix, in a suitable basis, is numerically a small perturbation from a diagonal 
real matrix whose nonzero entries are ±1. 

V Discussion 

We have observed, somewhat as a straw man argument, that at present it is numerically 
consistent to suppose that the neutrino mixing matrix is almost hermitian. 

If the neutrino matrix were hermitian, then it would be related by a change of basis to 
the matrix d = diag(-|-l, —1, —1). Moreover, despite having a large iKsl, a hermitian mixing 



matrix would necessarily conserve CP in oscillations [Eq. (III. 4)] 



Since tightening the upper bound on V^3 falsifies the hermitian ansatz, we have also presented 



a simple 1-parameter real parameterization of deviations from hermiticity [Eq. (IV.l)], which 



is useful for classifying perturbations away from the relation Ve3 = K-i- It is also easy to 



include the CP-violating PMNS angle in this parameterization [Eqs. (IV. 5) and (IV. 6)]. 



The lesson here is not whether the condition V"'' = V^ is approximately true, but rather 
that it is still experimentally consistent to perturb around an ansatz that is significantly 
different from the often-studied tribimaximal mixing matrix. We emphasize that one quan- 
titative measure of whether an ansatz is a good starting point is the 5'0(3)-invariant angle 
from a suitably extracted "snapshot" of the data [Eqs. (II. 8) and ( II.13[ )]. 



Finally, we invite the reader to find a possible theoretical origin for the condition V'' ~ V. 

Note added: After our work was submitted, W. Rodejohann called to our attention an 
earlier study [H] of the condition \V\'^ = \V\. This also brought to our attention a paper by 
Joshipura and Smirnov [12], which points out that the condition \V\'^ = \V\ can be obtained 
in models for which W Mu^/IA* = D^^d/ and U^M^U = D^, in a. self-evident notation, so 
that to leading order |Vckm| = I and |Vpmns| = |W^W|. 
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